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1 . 


SCOPE 




This document contains a design specification for implementing 
two Procedure 1 follow-on features. These features are; 

• On option, determine the best k of N passes in feature 
selection. 

• On option, include a-priori weighting based on pixel popula- 
tion in clusters in feature selection. 

These additions will be made to the SELECT processor of the EOD- 
LARSYS system. 

EOD-LARSYS is operational on the UNIVAC 1108 EXEC II computer 
system located in Building 12. The system is batch oriented, 
and operated and maintained according to IDSD procedures. This 
system is currently being converted to the Purdue-LARS 370/148 
system. 

This document assumes the reader is familiar with both Procedure 
1 and the EOD-LARSYS system. 




2. APPLICABLE DOCUMENTS 


• Final Design Specification for EOD-LARSYS Procedure 1, 
Houston, Texas, August 1977, JSC-12742, LEC-10417. 

• Job Order: 63-1347-1695 




3 . SYSTEM DESCRIPTION 



3.1 HARDWARE DESCRIPTION 
N/A 

3 .2 SOFTWARE DESCRIPTION 

Two new features have been added to the SELECT processor of the 
EOD-LARSYS system irk support of Procedure 1 follow-on require- 
ments . 

The first is the selection of the best k of N passes based on 
the overall separability criterion selected, where N is the 
total number of passes accounted for in the CHANNELS input and 
k is any number less than N. The number of channels per pass 
can be specified. 


The second is the modification of intersubclass weights by a 
multiplier based on the number of pixels in each subclass: 



where TOTSUB is the number of subclasses and N(i) is the number 
of pixels in subclass i. 


The following example illustrates the besk k of N selection, 
together with a-priori weighting. 

Control Cards : 


1, 2, 3, 4, 5 
UNIT=1, PILE=1 

1, 2, 4, 5, 6, 7, 9, 10, 11, 12, 14, 15 





SUBCLA 

STATFI 

CHANNEL 


PROCEED 6 
BSPASS 2 
APRIOR 

NCPASS 5 
*END* 


Here the number of channels per pass is 5, 3 passes are included 
(deduced from channels included ) , and the best 2 passes are to 
be determined. 

Each combination of 2 of 3 passes is considered. These are 


(1) 


2, 

4, 

5, 6, 7, 9, 10 

(2) 


2, 

4, 

5, 11, 12, 14, 15 

(3) 

6 , 

7, 

9, 

10, 11, 12, 14, 15 


For each, the overall separability measure is computed, and the 
optimum set of features (channels) is selected and output, 
through the INFORM labeled common block, this set is made 
available to the CLASSIFY processor. 

3.2.1 SOFTWARE COMPONENT NO. 1 (SELECT ) 

The SELECT subprogram is the main driver for the SELECT processor. 

3. 2. 1.1 Linkages 

Subprogram SELECT is invoked from the MONTOR monitor routine 
upon reading the $SELECT control card. It calls subprograms 
SETUP 4, ORDER, PRELIM, EXSRCH, WHRPLC , USERIN , GENRPT, PLOT, 
EVLFET. 

3 . 2 . 1 . 2 Interfaces 

SELECT interfaces with other routines through a calling sequence, 
and common blocks FSL, GLOBAL, INFORM, and BESTKN. BESTKN is 
newly incorporated for these options, it is coded as 


COMMON /BESTKN/ KPPPTS(60), IPRIOR, KBEST, NCPASS 


3.2 .1.3 Inputs 
Calling Sequence: 

Subr. SELECT (ARRAY, TOP) 


Parameter 

Dimension 

In/Out 

Description 

ARRAY 

TOP 

In/out 

Utility storage 


(presently set 
= 10500) 


for various arrays 

TOP 

1 

In 

Dimension of 




vector ARRAY 


3. 2. 1.4 Outputs 
N/A 

3. 2. 1.5 Storage 

Code; 1602g Data: 30337g 

3. 2. 1.6 Description 

The changes to subprogram SELECT involve computation of features 
(channels) to exhaust all combinations of k of N passes. Data 
vectors INDPER and PERM have been added to assist. Common block 
BESTKN has been added. For each combination of channels, sub- 
programs EVLFET and GENRPT are invoked for separability measure 
computation and results printout. 

3. 2. 1.7 Flowchart 
N/A 


3. 2. 1.8 Listing 
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3 -2 -2 SOFTWARE COMPONENT NO. 2 (SETUP4) 


Subprogram SETUr4 roads the control cards and initializes data 
and option switches for the SLECT processor. 

3. 2. 2.1 Linkages 

SETUP4 is called by routine SELECT and calls subprograms NUMBER, 
ORDER, NXTCHR, FIND, WGTSLN, CRDSTA, GRPSCN, BMFIL, REDSAV, 
BSTCHK, PRTFLD, and WGTCHK. 

3 » 2 . 2 . 2 Interfaces 

SETUP4 interfaces with other routines through a calling sequence 
and common blocks INFORM, GLOBAL, FSL, and BESTKN. 

3 . 2 . 2 . 3 Inputs 
Calling sequerce: 

Subr. SETUP4 (ARRAY, TOP, STOPFG, JTIME, SUBRAY, SUBSIZ) 


Parameter 

Dimension 

In/Out 

Description 

ARRAY 

TOP 

In/Out 

Utility storage for 
various arrays 

TOP 

1 

In 

Dimension of ARRAY 
vector 

STOPFG 

1 

Out 

Stop switch, set equal 
to 1 when $END* control 
board is read 

JTIME 

1 

Out 

Counter for number of 
times SELECT processor 
is invoked 


Parameter 

Dimension 

In/Out 

SUBRAY 

SUBSIZ 
(= 12000) 

Out 

SUBSIZ 

1 

In 

Now Control 

Cards ; 


APRIOR 

(Default- 
omit card) 


BSPASS 

N 

(no default) 


NCPASS 

N 



(default 4) 

Revised Control Cards: 

PROCED N 

(no default) 

SETUP4 inputs a statistics file 
reduction file (BMFIL) . 


Description 

Utility storage vc'' tor 

Dimension of vector 
SUBRAY 


This card sets the switch 
to modify inter-subclass 
weights 

N is the number of passes 
to be included in the 
best set 

N = number of channels 
per pass (acquisition) 


Set N = 6 to invoke the 
best k of N option 

(SAVTAP) (or cards) and a feature 


3. 2. 2. 4 Outputs 
. N/A 

3 . 2 . 2 . 5 Storage 

Code: 1772g Data: 1460g 

3 . 2 . 2 . 6 Description 

Subprogram has been modified to handle three additional control 
cards and an expanded option list for existing control card 




PHOCED and print these options as part of the ordinary user input 
summary. Control card data or switches are passed through 
variables IPRIOR, KBEST and NCPASS of common block BESTKN. 


'h. 1.2.1 Flowchart 
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3.2.3 SOFTWARE COMPONENT NO. 3 (PRELIM) 


Subprogram PRELIM sets up the separability measures for the full 
sot of channels specified by the user. This measure is either 
average divergence , transformed divergence , or Bhattacharriya 
distance. It also computes , an option , the inter-subclass 
weighting factors based on number of pixels per cluster and 
multiplies the weights by those factors. 

3.2. 3.1 Linkages 

Subprogram PRELIM is called by routine SELECT and calls DIVERG, 
AVEDIV, TRNDIV, BHTCHR, RWRITE, and ARITHMETIC routines SQRT 
and DEXP. 

3. 2. 3*2 Interfaces 

PRELIM interfaces with other routines through a calling sequence 
and common blocks IJFORM, FSL, and BESTKN. 

3 . 2 . 3 . 2 Inputs 
Calling Sequence: 

Subr. PRELIM (C^VMTX, AVEMTX, DIVTAB, WEIGHT, S, WRKRY, WRKSIZ) 


Parameter 

Dimension 

In/Out 

Description 

COVMTX 

(VARSZ2,NOCLS2) 

In/Out 

Subclass covariance 
matrices (lower 
triangular) 

AVEMTX 

(NOFET2 ,N0CLS2) 

In/Out 

Subclass mean vectors 

DIVTAB 

(DIVSTI2) 

In/Out 

Double Precision 
In ter- subclass 
divergences or dis- 
tances 

WEIGHT 

(DIVSI2) 

In/Out 

Inter-subclass 

weights 


T) 


Parameter 

Dimension 

In /Out 

Description 

S 

(VARS22,NOCLS2) 

Out 

"S'* matrices 




computed if CRIKEY 




* 1 

WRKRY 

(WRKSIS) 

In/Out 

Double Precision 


storage for 
feature subset 
covariance matrices 
and mean vectors 

WRKSIZ 1 In Computed as 12000 - 

SBASB in SELECT. 

SBASE « 1 if 
CRIKEY #1 * 2+ 
N0CLS2*VARSZ2 if 
CRIKEY*=1 


3 .2 . 3 .4 Outputs 
N/A 

3.2 .3.5 Storage 

Code: 671g Data; 1664 g 

3 . 2 . 3 . 6 De scription 

The changes to subprogram PRELIM involve the computation of the 
weighting factors and subsequent multiplication of the weights. 
These factors are computed from the KPPPTS array in newly-added 
labeled common /BESTKN/ which holds the number of pixels per 
cluster. 

3. 2. 3. 7 Flowchart 
N/A 


3. 2. 3. 8 Listing 
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3,2.4 SOFTWARE COMPONENT NO. 4 (vlTSTAT ) 


Subprogram GTSTAT acquires the reduced feature covariance matrices 
and moan vectors for each subclass. 

3.2. 4.1 Linkages 

Subprogram GTSTAT is called by EXSRCH, WHRPLC, FINT, USERIN, 
EVLFET, and calls subprogram TRNSFR. 

3 . 2 . 4 . 2 Interfaces 

GTSTAT interfaces with other routines through a calling sequence 
and common blocks INFORM and FSL. 

3.2. 4.3 Inputs 
Calling Sequence: 

Subr GTSTAT (COVMTX, AVEMTX, S, C0VMT2 , AVEMT2 , S2 , VEC, BMAT, 
WRKRY, IWRKSZ) 


Parameter 

Dimension 

In/Out 

Description 

COVMTX 

(VARSZ2,N0CLS2) 

In 

Covariance matrices 

AVEMTX 

(NOFET2 ,N0CLS2) 

In 

Mean vectors 

S 

(VARS Z 2 ,N0CLS2) 

In 

S matrices 

COVMT2 

(VARSZ4,NOCLS2) 

Out 

Double Precision 
reduced covariance 
matrices 

AVEMT2 

(N0FET4 ,N0CLS2) 

Out 

Double Precision 
reduced mean vectors 

S2 

(VARSZ4,NOCLS2) 

Out 

Double Precision 
reduced S matrices 

VEC 

N0FET4 

In 

Reduced channel set 


Parameter 


BMAT 


WRKRY 

IWRKSZ 


Dimension 
NOPET2 *NOFET4 


IWRKSZ 

1 


In/Ou t Description 

In Double Precision 

B-matrix for linear 
combination of 
channels (PRICKEY 
= 3 or 4) 

In/Out Double Precision 

In See WKRSIZ parameter 

PRELIM subprogram 


3 . 2 . 4 . 4 Outputs 
N/A 

3 . 2 . 4 . 5 Storage 


Code: 456 q Data: 104„ 

o 8 


3. 2. 4. 6 Desci'iption 


The change to GTSTAT consists of a transfer if PRCKEY=G to the 
channel reduction section. This is the same location reached 
for PRCKEY = 1,2, and 5. 


3.2.'’. 7 Flowchart 
N/A 


3. 2. 4. 8 Listing 
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3.2.5 SOFTWARE COMPONENT NO. 5 (GENRPT) 


Subprogram GENRPT handles the reporting of results through 
printer output. 

3. 2. 5.1 Linkages 

GENRPT is called by routine SELECT, and calls UNIVAC RANDIO 
routine RREAD. 

3 . 2 . 5 . 2 Interfaces 

GENRPT interfaces with other routines through a calling sequence 
and common blocks GLOBAL, FSL, and INFORM. 

3.2 .5. 3 Inputs 
Galling Sequences 

Subr. GENRPT (CLSNAM, WEIGHT, PIVTAB, WRKRY, IWRKSZ, FEJVEC) 


Parameter 

Dimension 

In/Out 

Description 

CLSNAM 

(NOCLSZ) 

In 

Class names 

WEIGHT 

(DIVISZ) 

In 

inter-subclass 

weights 

DIVTAB 

(DIVSIZ) 

In 

Double Precision 

inter-subclass 

separabilities 

WRKRY 

IWRKSZ 

In 

Double Precision 
separabilities for 
full feature set 

IWRKSZ 

1 

In 

See WRKSIZ parameter 
of routine PRELIM 

FETVEC 

30 

In 

Selected features 



3. 2 .5. 4 Outputs 

Printed reports are outputted. 


3 . 2 . 5 . 5 Storage 

Code: 1116g Data: 523g 

3 • 2 . 5 . 6 Description 

The changes to GENP.PT consist of privisions for extra printout 
if PRCKEY=6. 

3. 2. 5. 7 Flowchart 
N/A 

3. 2. 5. 8 Listings 
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3.2.6 SOFTWARE COMPONENT NO. 6 (REDDAT) 


Subprogram REDDAT is one of a set of utility routines involved 
with reading statistics files. 

3 . 2 . 6 * 1 Linkages 

REDDAT is called by REDSAV. 

3. 2. 6. 2 Interfaces 

REDDAT interfaces with other routines through a calling sequence 
and common blocks INFORM, GLOBAL and BESTKN. 

3 . 2 . 6 . 3 Inputs 
Calling Sequence: 

Subr. REDDAT (COVAR, AVAR, CLSDES , SUBNO, SUBDES , FLDSAV, VERTEX, 
COV, AVEN, CLSDS, SUBNOS, SUBDS, FLDSV, VEPTX, 
NOFEAT, VARSIZ, NOCLS, NOFLD, NOSVB, FETVEC) 

These parameters are all associated with statistics file (SAVTAP) 
information . 

3 . 2 . 6 . 4 Outputs 
N/A 

3 . 2 . 6 . 5 Storage 

Code: 1014g Data: 212g 

3. 2. 6. 6 Description 

The additions to REDDAT involve the storage of the number of 
points per subclass (after grouping) into array KPPPTS of the 
newly-added common block BESTKN. 



3 . 2 . 6 c 7 Flowchart 
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4 . OPERATION 


For the duration of the UNIVAC 1.108 EXEC 2 syctem implementation, 
these options are included in the EOD-LARSYS system by means of 
PCF tapes X23778 and X14135. One of these tapes should be 
referenced as unit Z in run setups. 

These options are also included in the EOD-LARSYS implementation 
on the Purdue-LARS 370/148 system. Detailed instructions on the 
usage of this system will be provided at a later date. 




